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Abstract 

, — . I investigate the elastic electron scattering off nuclei far from the stability line. The effects of the 

f— ^ . neutron and proton skins and halos on the differential cross sections are explored. Examples are 
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given for the charge distribution in Sn isotopes and its relation to the neutron skin. The neutron 



(^ ' halo in ^^Li and the proton halo in ^B are also investigated. Particular interest is paid to the 

fvq ■ inverse scattering problem and its dependence on the experimental precision. These studies are of 

particular interest for the upcoming electron ion colliders at the GSI and RIKEN facilities. 
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I. INTRODUCTION 

The use of radioactive nuclear beams produced by fragmentation in high-energy heavy- 
ion reactions lead to the discovery of halo nuclei, such as ^^Li, about 20 years ago [l|. 
Nowadays a huge number of /3-unstable nuclei far from stability are being studied thanks 
to further technical improvements. Unstable nuclei far from stability are known to play an 
important role in nucleosynthesis. Detailed studies of the structure and their reactions will 
have unprecedented impact on astrophysics [2(]. 

The first experiments with unstable nuclear beams aimed at determining nuclear sizes 
by measuring the interaction cross section in high energy collisions [1]. Successive use of 
this technique has yielded nuclear size data over a wide range of isotopes. Other techniques, 
e.g. isotope-shift measurements, have allowed to extract the charge size. The growth of a 
neutron skin with the neutron number in several isotopes have been deduced from nuclear- 
and charge-size data |3|. 

The charge (nucleon) density distribution can also be determined by electron (hadron) 
scattering experiments. Among hadron scattering, proton elastic scattering at intermediate 
energies is a good tool to probe the nucleon density distributions, due to its larger mean-free 
path in the nuclear medium. But, undoubtedly, it is the electron scattering off nuclei that 
provides the most direct information about charge distribution, which is closely related to 
the spatial distribution of protons [4]. 

A technical proposal for an electron-heavy-ion collider has been incorporated in the 
GSI/Germany physics program |5|. A similar program exists for the RIKEN/ Japan fa- 
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In both cases the main 



cility using Self-Confining Radioactive Ion Targets (SCRITs) 
purpose is to study the structure of nuclei far from the stability line. The advantages of using 
electrons in the investigation of the nuclear structure are mainly related to the fact that the 
electron-nucleus interaction is relatively weak. For this reason multiple scattering effects are 
usually small and the scattering process is described in terms of perturbation theory. Since 
the reaction mechanism in perturbation theory is well under control the connection between 
the cross section and quantities such as charge distributions, transition densities, response 
functions etc., is well understood [Tj]. 

Theoretically, the shape of the density distribution includes detailed information on the 
internal nuclear structure. In the independent particle shell model the density distribution 



is the squared sum of the single-particle wave functions. No measurement of either nucleon 
distribution or the charge distribution for short-lived radioactive nuclei has been made so 
far. 

Under the impulse approximation, or plane wave Born approximation, the charge form 
factor can be determined from the differential cross section of elastic electron scattering. 
Since the charge distribution, Pch{f), is obtained from the charge form factor by a Fourier 
transformation, one can experimentally determine pch{f) by differential cross-section mea- 
surements covering a wide range of momentum transfer q. This leads to information on 
the size and diffuseness when the charge form factor is measured at least up to the first 
maximum. To accomplish this with a reasonable measuring time of one week, a luminosity 
larger than 10^^ cm~^s~^ is required, for example for the ^^^Sn isotope |5|. 

On the theoretical side the difference between the proton and neutron distributions can 
be obtained in the framework of Hartree-Fock (HF)method (see for example [8|) or Hartree- 
Fock-Bogoliubov (HFB) method (see for example J9|, Il0|). As a rule of thumb, a theoretical 
calculation of the nuclear density is considered good when it reproduces the data on elastic 
electron scattering. But some details of the theoretical densities might not be accessible 
in the experiments, due to poor resolution or limited experimental reach of the momentum 
transfer a. Recent works have also looked at electron scattering from halo nuclei, see e.g. 

iii,yu- 



refs. 



In this work I study the general features of elastic electron scattering off unstable nuclei. 
Here I focus on using the general features of what is theoretically known about skins and 
halos to look for their respective signatures in the data. Only a few specific and represen- 
tative nuclear density cases are used. The paper is organized as follows. After a general 
introduction of elastic electron scattering in section 2, I show in section 3 how the telltales of 
neutron and proton skins and halos become visible in the differential cross sections for elas- 
tic electron scattering. Section 4 presents a study of the inverse scattering problem, which 
poses a challenge for extracting the charge density profiles from the experimental data. The 
conclusions are presented in section 5. 



II. ELASTIC ELECTRON SCATTERING 

In the plane wave Born approximation (PWBA) the relation between the charge density 
and the cross section is given by 

d^ J PWBA 1 + (^E/Ma) sin^ (9/2) I ^'^ ^^^ I ' ^ ^ 

where (Jm = {e'^/4:E^) cos^ (9/2) sin^"^ (9/2) is the Mott cross section, the term in the denom- 
inator is a recoil correction, E is the electron total energy. Ma is the mass of the nucleus 
and 6 is the scattering angle. 

The charge form factor Fch {q) for a spherical mass distribution is given by 

POD 

Fch (q) =4:71 dr r^jo (qr) pch (r) , (2) 

Jo 

where q = 2k sin (6/2) is the momentum transfer, hk is the electron momentum, and E = 

^yh'^k'^c^ + mlc^. The low momentum expansion of eq. [2] yields the leading terms 

FcH{q)lZ = l-'^^{rl,) + .... (3) 

Thus, a measurement at low momentum transfer yields a direct assessment of the mean 
square radius of the charge distribution, (r^^) . However, as more details of the charge 
distribution is probed more terms of this series are needed and, for a precise description of 
it, the form factor dependence for large momenta q is needed. 

A theoretical calculation of the charge density entering eq. [2] can be obtained in many 
ways. Let pp (r) and p„ (r) denote the point distributions of the protons and the neutrons, 
respectively, as calculated, e.g. from single-particle wavefunctions obtained from an average 
one-body potential well, the latter in general being different for protons and neutrons. If 
fEp (r) and fEn (r) are the spatial charge distributions of the proton and the neutron in the 
non-relativistic limit, the charge distribution of the nucleus is given by 

PcH (r) = j Pp (r') JE, (r - r') d?r' + j p^ (r') Jeu (r - r') rfV. (4) 

The second term on the right-hand side of eq. H] plays an important role in the inter- 
pretation of the charge distribution of some nuclear isotopes. For example, the half-density 
charge radius increases 2% from ^°Ca to ^^Ca, whereas the surface thickness decreases by 
10% with the result that there is more charge in the surface region of '^''Ca than of ^^Ca 



|l4| . This also implies that the rins charge radius of ^^Ca is slightly smaller than that of 
^°Ca. The reason for this anomaly is that the added j-jj^ neutrons contribute negatively to 
the charge distribution in the surface and more than compensate for the increase in the rms 
radius of the proton distribution. 

For the proton the charge density ]ep if) i^ ^Q- SI is taken as an exponential function, 
corresponding to a form factor JE-p (q) = (1 + g^/A^)^^ (see Appendix A). For the neutron 
a good parametrization is fEnil) = —fJ'nTfEp{(l)/{^ +Pt), where /i„ is the neutron magnetic 
dipole moment and r = q'^/4mi^. We will use A^ = 0.71 fm~^ (corresponding to a proton 
rms radius of 0.87 fm) and p = 5.6, which Galster et al. [l5| have shown to reproduce 
electron-nucleon scattering data. 

Eqs. [IHare based on the first Born approximation. They give good results for light nuclei 
(e.g. ^^C) and high-energy electrons. For large-Z nuclei the agreement with experiment is 



only of a qualitative nature. The effects o 



distortion of the electron waves have been 



studied by many authors (see, e.g. ref. la, ll7|, llSl]). More important than the change in the 
normalization of the cross section is the displacement of the minima. It is well known that 
a simple modification can be included in the PWBA equation reproducing the shift of the 
minima to lower g's. One replaces the momentum transfer q in the form factor of eq. [T]with 
the effective momentum transfer ge// = Q' (1 + 3Ze^/2i?c/i-E'), where E is the electron energy 
and Rch — 1.2 A^^^ fm. This is because a measurement at momentum transfer q probes in 
fact |F(g)p at g = q^fj due to the attraction the electrons feel by the positive charge of 
the nucleus. This expression for q^jf assumes a homogeneous distribution of charge within 
a sphere of radius Rch- 

A realistic description of elastic electron scattering cross sections requires full solution of 
the Dirac equation. The Dirac equation for elastic scattering from a charge distribution can 



be found in standard textbooks, e.g. [19|. Numerous DWBA codes based on Dirac distorted 
waves have been developed and are public. Since the inclusion of Coulomb distortion is 
straight-forward, we will concentrate on the information which can be extracted from elastic 
electron scattering, which is encoded in the form-factor of eq. El Later we will study the 
problem of extracting the full nuclear density from the form factor. For simplicity, we 
will also neglect the effect of the charge distribution of the neutron itself. This effect is 
particularly important for large g-transfers. For heavy nuclei it can change \Fch{q)\ by 
10-20% in the q-range of 1.5 - 2 fm~^. 



III. SKINS AND HALOS 



Neutron Skins 



Appreciable differences between neutron and proton radii are expected [20| to characterize 
nuclei at the border of the stability line. The liquid drop formula expresses the binding 
energy of a nucleus with A^ neutrons and Z protons as a sum of bulk, surface, symmetry 
and Coulomb energies EjA = ay A - asA^/^ - S{N - Zf/A - acZ'^ j A^l'^ =f OpA^^/^, where 
ay, asi Op, 5* and ac are parameters fitted to experimental data of binding energy of nuclei. 
This equation does not distinguish between surface (S*) and volume {V) symmetry energies. 
As shown in ref. J2l|], this can be achieved by partitioning the particle asymmetry as 
N ~ Z = Ns ~ Zs + Ny — Zy. The total symmetry energy S then takes on the form 
S = Sv{Ny - ZyfjA + SsiJ^R- ZsfjA^I^. Minimizing under fixed N - Z leads to an 
improved liquid drop formula 2l(| with the term S{N — Z^'^ j A replaced by 



The same approach also yields a relation between the neutron skin /?„„ = i?„ — /?„, and 



5*5, Sy^ namely [21 1 



Rn-Rp _ A 7\- ^ N-Z-{ac/l2Sy)ZA^'^ 

R QNZ^ ^ ^' 6NZ l + {Ss/Sy)A^/^ ' ^^ 

where R = {Rn + Rp)/2, and Rn (Rp) is the half-density neutron (proton) distribution radius. 
Here the Coulomb contribution is essential; e.g. for N = Z the neutron skin Rnp is 
negative due to the Coulomb repulsion of the protons. A wide variation of values of 5*^ and 
5*5 can be found in the literature. These values have been obtained by comparing the above 
predictions for energy and neutron skin to theoretical calculations of nuclear densities and 



experimental data on other observables 
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22|. The values of ac = 0.69 MeV and Sy = 28 



MeV, which will be used here, are compatible with fits of the binding energy and symmetry 
energy of known nuclei J22|. Using these values and R = 1.2^1^/'^ fm, one gets, for large A 
(with A > 1, AZ ^ AV4), 



Rnp = Rn-Rp^ 0.8 ( -^ j {6 - 2.05 X IQ-^ZA-^/^) fm. 



(7) 



where S = {N — Z) j A is the asymmetry parameter. The ratio Sy/Ss varies within the range 
2 ^ Ss/Sy > 3 by adjusting eq. [6]to the dependencies of skin sizes on asymmetry and on the 



neutron-proton separation-energy difference constraints [22] . If one assumes that the central 
densities for neutrons and protons are roughly the same and that they are both described 
by a uniform distribution with sharp-cutoff radii, i?„ and Rp, one finds Rnp — O.SA^^^S fm. 
This shows that the sharp sphere model is too simple and that in this model the central 
densities for proton and neutrons have to differ if eq. is to remain valid. 



On the experimental front, a study of antiprotonic atoms published in reference [23 1 
obtained the following fitted formula for the neutron skin of stable nuclei in terms of the 
root mean square (rms) radii of protons and neutrons 

Ar„p = {rlY^^ - {rlY^^ = (-0.04 ± 0.03) + (1.01 ± 0.15) 5 fm. (8) 

The relation of the mean square radii with the half-density radii for a Fermi distribution is 
roughly given by (r^) = 3-R^/5 -|- Tvr^a^/S, where a„ is the diffuseness parameter. For heavy 
nuclei, assuming an = ap <^ Rn, Rp, one gets approximately the same linear dependence on 
the asymmetry parameter as in eq. [7] (neglecting the small contribution of the A~^^^ term). 

In this work we will use eq. [8] as the starting point for accessing the dependence of 
electron scattering on the neutron skin of heavy nuclei. Applying it to calcium isotopes as 
an example, one obtains that the neutron skin varies from —0.15 fm for '^^Ca (proton-rich 
with negative neutron skin) to 0.25 fm for ^^Ca. A negative neutron skin obviously means 
an excess of protons at the surface. For unstable nuclei no study as the one leading to eq. 
[8] is experimentally available. Nucleon knockout reactions with secondary unstable nuclear 
beams, and other techniques, have been used to determine the interaction radii of neutron 
rich nuclei Jj. These radii are thought to represent the extension of the neutron distribution. 
But very little is known about the charge distribution in neutron-rich nuclei. 

For heavy nuclei the charge and neutron distributions can be described by a Fermi dis- 
tribution. The diffuseness is usually much smaller than the half-density radius, Op „ -C Rp,n- 



The neutron skin is then given by Rnp = Rn — Rp — a/5/3 Ar„p. We assume further that the 
nuclear charge radius is represented by Rp = 1.2^4^/^ fm — Rnp/2 with Ar^p given by eq. [3 
Figure [1] shows the charge form factor squared for elastic electron scattering off tin isotopes, 
as a function of the momentum transfer. The two solid curves are for the extreme values 
of the asymmetry parameter 6 = {N — Z)/A, that is 5 = (A^ = Z = 50), and 6 = 4/9 
(A^ = 90). The curves form an envelope around other curves with intermediary values of S. 
The first and second minima of the form factors occur at qi = AA9/Rp and g2 = 7.73/ Rp, 
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FIG. 1: The charge form factor squared for elastic electron scattering off tin isotopes, as a function 
of the momentum transfer. The two curves are for the extreme values of the asymmetry parameter 
6 = {N- Z)/A, that is 6 = {N = Z = 50), and 6 = 4/9 {N = 90). The curves form an envelope 
around other curves with intermediary values of 6. 

respectively, corresponding to the zeroes of the transcendental equation tan {qRp) = qRp. 

Using the approximations discussed in the above paragraph, i.e. Rp = 1.2^4^/^ fm — i?„p/2 
with Avnp, Rnp = Afnp, and the experimental value given in eq. [HI we conclude that the 
linear dependence of Rp with the neutron skin (and with the asymmetry parameter 5) , also 
implies a linear dependence of the position of the minima, 
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92 = 1.72 gi . 



(9) 



For ^''^Sn the first minimum is expected to occur at gi = 0.806 fni = 159 MeV/c, while 
for ^32Sn it occurs at gi = 0.754 fm"^ = 149 MeV/c. 

Figure [2] shows the value of gi for calcium, tin and uranium isotopes, as a function of the 
asymmetry parameter 5. The variation of gi with the neutron skin of neighboring isotopes, 
Agi ~ 2/A^/^ fm~^, is too small to be measured accurately. The first minimum, gi, changes 
from 220 MeV/c for ^^Ca to 204 MeV/c for ^^Ca, approximately 7%, which is certainly 
within the experimental resolution. Of course, sudden changes of the neutron skin with 5 
might happen due to shell closures, pairing, and other microscopic effects. 
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FIG. 2: Position of the first minimum qi for elastic electron scattering off calcium, tin and uranium 
isotopes, as a function of the asymmetry parameter 5. 

To be more specific, let us assume that a reasonable goal is to obtain accurate results for 
the charge radius (r^) , so that S (r^) < 0.05 fm. This implies that the measurement of 
gi has to be such that (Agi/gi) < gi [fm~ ] %, with gi in units of fm~^ and the right-hand 
side of the inequality yielding the percent value. For ^^Ca, one has qi = 1.11 fm~^, meaning 
that the experimental resolution on the value of gi has to be within 1% if 6 (rT) < 0.05 
fm is a required precision. The situation improves for heavier nuclei, as becomes evident 
from figure O 

The neutron skin dependence is also seen in the height of the second bump, after the 
first minimum. This bump occurs at g^ — Q/Rp. For a Woods-saxon density this peak 
will be reduced compared to the first maximum by a factor [g{qma)/ 12]"^, where g{qmCi) is 
a function of the diffuseness parameter a. g is closely given by the Fourier transform of 
an Yukawa function, its value at Qma is of the order of 1/2 and its variation around g^ is 
weak, g{qa) ~ 1/ (1 + g^a^) (see Appendix A). Thus, the dependence of the height of the 
second maximum upon the neutron skin is a less appropriate tool than the location of the 
first minimum. Of course, the ultimate test of a given theoretical model will be a good 
reproduction of the measured data, below and beyond the first minimum. 



B. Neutron halos 

For light halo nuclei composed of a core nucleus and an extended distribution of halo 
nucleons, the nuclear matter form factor can be fitted with the simple expression 

F{q)/A ={l-g) exp {-q^l/A) + — ^, (10) 

where g is the fraction of nucleons in the halo. In this expression the first term follows from 
the assumption that the core is described by a Gaussian and the second term assumes that 
the halo nucleons are described by an Yukawa distribution (see Appendix A). Taking ^^Li as 
an example, the following set of parameters can be used g = 0.18, ai = 2.0 fm and 02 = 6.5 
fm. This means that g = 2/11 nucleons are in the halo, the size of the core is roughly 2 fm, 
and the size of the halo is 6.5 fm. Although only few nucleons are in the halo they change 
dramatically the shape of the form factor, as shown in figure [3l The dashed-dotted curve is 
the squared charge form factor. The dotted curve applies for the matter density of the halo 
neutrons. When the core (dashed curve) and halo nucleon distributions are combined the 
squared form factor for the total matter distribution (solid curve) clearly displays the halo 
signature. Thus, even when the individual contribution of the halo nucleons is small and 
barely visible in a linear plot of the matter distribution, it is very important for the form 
factor of the total matter distribution. It is responsible for the narrow peak which develops 
at low momentum transfers. This signature of the halo was indeed the motivation for the 
early experiments with radioactive beams. The narrow peak was observed in momentum 
distributions following knockout reactions ij]. 

Elastic electron scattering will not be sensitive to the narrow peak of |-F(q')| (matter 
distribution form factor) at small momentum q, but to the the charge distribution form 
factor, |Fc/i(g)| , which in the case of ^^Li will be similar to the dashed-dotted curve in 
figure [3l The determination of this form factor will tell us if the core has been appreciably 
modified due to the presence of the halo nucleons. It is worthwhile mentioning that many 
of what we call core nuclei are also short-lived and very little is known about their charge 
distribution. 

In order to explain the spin, parities, separation energies and size of exotic nuclei con- 
sistently a microscopic calculation is needed. One possibility is to resort to a Hartree-Fock 
(HF) calculation. Unfortunately, the HF theory cannot provide the predictions for the sep- 
aration energies within the required accuracy of hundred keV. Here we use a simple and 
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FIG. 3: Form factor, \F{q)\'^ /A^, for ^^Li showing the individual contributions of the core and of 
the halo nucleons. The dotted (halo), dashed (core) and solid (total) curves are form factors for 
the matter distribution. The dashed-dotted curve is the charge form factor, \Fch{q)\'^ /Z'^- 



tractable HP method 25| to generate synthetic data for the charEre-distribution of ^^Li. 
Detai. of th. method . dLcnbed . .. U. A.u.i„. .pheHcal =.„„... t.e e,.at.n 

for the Skyrme interaction can be written as 



-V- 



-V + V{r) 



^pi{r) = ei^pi{r) 



(11) 



2m* [r) 

where m*{r) is the effective mass. The potential V{r) has a central, a spin-orbit, and a 
Coulomb term 

I' (r) = Vcentral + Vspin-orbit + ^Coulomb ■ (12) 

The central potential is multiplied by a constant factor / only for the last neutron config- 
uration: 

/ 7^ 1 for last neutron configuration 

/ = 1 otherwise. 

Thus, the last neutron configuration (last orbits) is treated differently from the other orbits 
in the HF potential in order to reproduce the neutron separation energy of the neutron-rich 
nucleus. The factor / in Eq. [13] is arbitrary (/ = 0.82 for the last neutrons in ^^Li [25]). It 
roughly scales with the inverse of the fraction of nucleons in the halo, and simulates a weaker 



Vr 



central 



(r) = fVHFir), 



(13) 
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potential at the halo region. This model was successful to explain most features of the light- 



neutron rich nuclei 



25 



271] . It can also explain the magnitude of the nuclear sizes. In order 
to obtain the nuclear sizes, the rms radii of the occupied nucleon orbits are multiplied by 
the shell model occupation probabilities, which are also obtained in the calculations. The 
final radius is obtained by adding the core radius. It is important to notice that the physics 
of ^^Li is not treated very well because of the pairing interaction. This is needed to make 
^^Li bound while ^°Li is unbound. In this aspect, the model adopted here is only useful as 
a qualitative tool. 



Jhalo yFIcal (fm) x/FIexp i^"^) 



9Li (core) 2.45 2.43 ± O.OTt 



11 



Li 



lPl/2 


5.36 


2Si/2 


7.61 




3.26 



SM 3.26 3.62 ± 0.19+ 

^ from ref. [29] . 

Table 1 - Single particle properties of Li. The second column gives the spins of the most probable 
occupied orbits. The third column is the result of HF calculations for the rms radii associated with 
these orbits, and the last column gives the rms radii of the matter distribution of these nuclei. 

As the effective interaction, a parameter set of the density dependent Skyrme force, so 
called BKN interaction [30|, is adopted. The parameter set of BKN interaction has the 
effective mass m*/m =1 and gives realistic single particle energies near the Fermi surface in 
light nuclei. The original BKN force has no spin-orbit interaction. In the present calcula- 
tions, I introduce the spin-orbit term in the interaction so that the single-particle energy of 
the last neutron orbit becomes close to the experimental separation energy. In this way, the 
asymptotic form of the loosely-bound wave function becomes realistic in the neutron-rich 
nucleus. The large r.m.s. radii of the valence neutron orbits is attributed to the small sep- 
aration energy. The calculated value is enough to create the halo structure in the HF wave 
functions. In ^^Li, the last occupied orbits is taken to be lpi/2 and 2si/2. 

We took the cut-off radius of H-F calculation to be i? = 40 fm which is necessary to 
include properly the loosely bound nature of the neutron wave function. The separation 
energy for the valence nucleons is S2n{theo) = 0.274 MeV, which is to be compared with 
the experimental one S2n(ea;p) = 0.247 ±0.080 MeV 



28|. The column indicated by jhaio in 
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Table 1 displays the most probably occupied orbits. The final radius is obtained by adding 
the core radius, and is given in the row indicated by SM. 

The elastic form factor for the matter distribution obtained in the HF calculations are 
very close to the ones calculated by the empirical formula (TDl In figure [3] we show the 
charge form factor, |Fc/j(g)| , by a dashed-dotted line. There is very little difference from 
the simpler empirical parametrization of eq. [101 The lack of minima, and of secondary 
peaks (as in eq. [TOl) . makes it difficult to extract from \Fch{q)\ more detailed information 
on the charge-density profile. Parametrizations like eq. [10] which apparently lack of a sound 
physical basis are not particular to loosely bound nuclear systems. For example, in the case 



of ^Li a good fit to experimental data was obtained with 



3l| 



\Fch{q)\ oc exp(-a q ) - Cq exp(-6 q 



with a = 0.933 fm, b = 1.3 fm, and C = 0.205. However, the data 311] cannot be fitted by 



using a model in which the nucleons move in a single-particle potential. 

C. Proton Halos 

Here I will consider ®B as a prototype of proton halo nucleus. This nucleus is perhaps 
the most likely candidate for having a proton halo structure, as its last proton has a binding 
energy of only 137 keV. The charge density for this nucleus can be calculated in the frame- 



work of the Skyrme HF model. We will use here the results obtained in ref. 32J], where 



axially symmetric HF equations were used with SLy4 [33|| Skyrme interaction which has 
been constructed by fitting the experimental data on radii and binding energies of symmet- 
ric and neutron-rich nuclei. Pairing correlations among nucleons have been treated within 
the BCS pairing method. The form factor squared for the charge density in ^B is shown 
in figure HJ normalized to mass number for matter distribution, and to charge number for 
charge distribution. 

The width of the charge form factor squared corresponding to the plot in figure [His Ach = 
0.505 fm^^ = 99.6 MeV/c. The corresponding values for the form factors for neutron and 
total matter distributions are, respectively, A„ = 0.512 fm^"*^ = 101 MeV/c and A^t = 0.545 
fm~^ = 108 MeV/c. This amounts to an approximately 10% difference between matter and 
charge form factors in ^B. 
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FIG. 4: Form factor squared for ^B showing the contribution of neutron, \F{q)\'^ /N'^, (dashed) 
charge, \F{q)\'^ /Z"^, (dotted) and total, |F(g)p/j4^, (solid) matter distribution of nucleons. 

The proton halo in ^B is mainly due to the unpaired proton in the P3/2 orbit. It is clear 
that for such a narrow halo the size of the nucleon also matters. The relevance of the nucleon 
size is shown in figure |5] where a slice of the nucleon is included in a thin spherical shell of 
radius r and thickness dr from the center of the nucleus. If the position of the nucleon is 
given by R, the part of the proton charge included in the spherical shell is given by 

dpch = '^T^r'^dr / d9 Pp(x)sm9, (14) 

Jo 

where pp (x) is the charge distribution inside a proton at a distance x from its center. The 

coordinates are shown in figure [5l They are related by x'^ = r'^ + R"^ — 2rRcos6. The 

contribution to the nuclear charge distribution from a single-proton in this spherical shell is 

thus given by Afp {R, r) = dpch/4:Trr^dr. 

Assuming that the charge distribution of the proton is described either by a Gaussian or 

a Yukawa form, the integral in eq. [T3] can be performed analytically, yielding 



AfP{R,r) 



1 
47rV2ari? 



exp 



(R-rf 



exp 



\R-r\ 



— exp 

— exp 



(R+rf 



a2 

l-R+rl 



for a Gaussian dist. 

for a Yukawa dist., 

(15) 

where a is the proton radius parameter. 

The charge distribution at the surface of a heavy proton- rich nucleus, 6pch (r), may be 
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described as a pile-up of protons forming a skin. Let Ui be the number of protons in the 
skin and Ri their distance to the center of the nucleus. One gets 

5pch{r) = ^niMp{r,Ri). (16) 

i 

Assuming Ri to be constant, equal to the nuclear charge radius i?, and using eqs. [TB] it 

is evident that while the density at the surface increases, its size R and width a, remain 

unaltered. The form factor associated with this charge distribution is given by 

^^ ,. 47r ^-^ f°° , A/-/ 7-.\ • / \ / 2\ v^ sin(g_Rj) 

oF^h [q) = — y ,^i I dr r Np (r, Ri) sm [qr) = exp [-qa ) } rii , (17) 

q ^ Jo ^ ' ^ qR^ 

where the last result is for the Gaussian distribution. An analytical expression can also be 

obtained for the Yukawa distribution. For Ri = R, expression [17] shows that the increase 

of density in the skin does not change the shape of the form factor, or of the cross section, 

but just its normalization. The decrease of the form factor with q is determined by a alone, 

and not by rii. If the charge of additional protons is distributed homogeneously across the 

nucleus including the skin, the form factor will not change appreciably, except for a small 

change in R. 

For a proton halo nucleus it is more appropriate to replace Yli ^i ~^ 

4:71 f dR R^ Afp (r, R) 6pch (R), where 6pch (-R) is the density change created by the extended 

wavefunction of the halo protons. We then recast eq. [T7] in the form 

5F,h (q) = — exp {-qa'') / dR R Sp^h (R) sin (qR) . (18) 

Q Jo 

The shape of the form factor is here dependent not only on the proton size but also on the 
details of the halo density distribution. For ^B, the halo size is determined by the valence 
proton in a P3/2 orbit. The density 6pch (R) due this proton can be calculated with a Woods- 
Saxon model. Using the same potential parameters as in ref. (34| we show in fig. [6] the form 
factor \6Fch {q)\ compared to the charge form factor \Fch {q)\ of figure HI It is evident that 
the halo contributes to a narrow form factor. However, in contrast to the neutron halo case 
shown in figure [3l the charge form factor of ^B does not show a pronounced influence of the 
halo charge distribution. 

The low energy expansion of the form factor, eq. [3|, allows the extraction of the rms 
radius of the charge distribution from 

dF,f, 



(d> = -6 



dig') 



(19) 
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FIG. 5: A spherical shell with radius r from the center of the nucleus will have a contribution from 
the charge inside a nucleon located at a distance R. 

Applying this relationship to the charge form factor used in figure E] we get (r^^^) = 2.82 
fm which is close to the experimental value {t1i^)J = 2.82 ± 0.06 fm [35]. The shape of the 
charge form factor can also be described by a Gaussian distribution with radius parameter 
a = 2.30 fm. In contrast to the case of ^^Li seen in figure [3l the proton halo in ^B does 
not seem to build up a two-Gaussian shaped form factor. This observation also seems to be 
in contradiction with the momentum distributions of ^Be fragments in knockout reactions 
using ^B projectiles in high energy collisions 2J]. Electron scattering experiments will help 
to further elucidate this property of proton halos. 



D. Comparison with inelastic scattering 

Loosely bound nuclei easily undergo breakup under any excitation. Such inelastic pro- 
cesses will always accompany elastic scattering. We follow the equations obtained in ref. 



36| to make a crude estimate of the effects of inelastic electron scattering off loosely bound 



nuclei. 

The inelastic cross section for the electro-excitation of loosely bound nuclei is given by 



36| 



da\ 



incl 



48^2' 



--eff 



1 2 



p 1 Vs~AE,-Sn 



,3/2 



dQdE.. 



n K^fic^ q^ 



Et 



(20) 



where E^ is the excitation energy, Vt denotes the scattering solid angle, e is the electron 
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a^ 0.5 - 







0.0 0.1 0.2 0.3 0.4 0.5 



FIG. 6: The charge form factor squared of ^B (sohd hne), |-Fj!^*"'(g)p/5^. Also shown is the charge 
form factor squared due only to the valence halo proton, |F^^ (g)p/l^. 

charge, e^g is the electric effective charge for the nuclear response (I will assume for simplicity 
that e[J=Ze), and Sn is the nuclear breakup energy, p is the electron momentum, q is the 
momentum transfer, and fi is the reduced mass, where it is assumed that the nucleus is 
composed of two clusters. 

We now integrate this equation over the excitation energy E^, and obtain 



da-. 



incl 



3^2 Z 



2„4 



(21) 



dn 4 SnHc^sm\e/2)' 

This equation is valid for large electron energies, so that E ^ meC^, and for small scattering 
angles, for which the approximation hq ~ 2psin(^/2) can be used. In fact, the minimum 
momentum transfer for an excitation energy E^ is given by gmin = ^k = E^/hc. Thus, the 
inelastic cross section tends to a maximum value at very small scattering angles (in contrast, 
the elastic scattering cross section increases indefinitely at small angles). Equation [2T] gives 
the behavior of the inelastic cross section when it starts to decrease with angle, shortly after 
the maximum at zero degrees. A characteristic feature emerging from eq. [21] is that the 
inelastic cross section is proportional to the inverse of the separation energy Sn- 

From equation [H neglecting nuclear recoil, one gets for the elastic scattering cross section 
dae\/dVL = Z^e'^/AE^ sm'^{9/2). The ratio between elastic and inelastic cross sections for 
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small scattering angles is thus given by 

dad 1 Snfic^ 1 



(22) 



This shows that the elastic scattering dominates over inelastic scattering for 6 <^ Omax ~ 
y^Snfic'^/E. Adopting typical values, i.e. E = 100 MeV, S'„ = 1 MeV and /ic^ = 10^ MeV, 
one gets Omax ~ 1/3 radians. For E = 1 GeV this value reduces to ^max ~ 1/30 radians. 
These are kinematical constraints which have to be taken into account in future electron 
scattering experiments. 

IV. INVERSE SCATTERING PROBLEM 

In PWBA the inverse scattering problem can be easily solved. It is possible to extract 
the form factor from the cross section and then, with an inversion of the Fourier transform, 
to get the charge density distribution 

1 f°° 

Pch (^) = TT^ / dq q^io (qr) F^h (g) • (23) 

^^ Jo 

As we discussed in previous sections, the PWBA approximation can be justified only for 

light nuclei (e.g. ^^C) in the region far from the diffraction zeros. For higher Z values the 

agreement with experiment is only of a qualitative nature. 

It is very common in the literature to use a theoretical model for p^h (^), e.g. the HF 
calculations discussed in the previous sections and compare the calculated F{q) with ex- 
perimental data. When the fit is "reasonable" (usually guided by the eye) the model is 
considered a good one. However, whereas the theoretical pch (^) can contain useful informa- 
tion about the central part of the density (e.g. bubble-like nuclei, with a depressed central 
density), an excellent fit to the available experimental data does not necessarily mean that 
the data is sensitive to those details. The obvious reason is that short distances are probed 
by larger values of q. Experimental data from electron-ion colliders will suffer from limited 
accuracy at large values of g, possibly beyond g = 1 fm^^. Thus it is useful to identify what 
are the conditions for reproducing the nuclear density within a theory independent fit. 

In order to obtain an unbiased "experimental" pch (^) one usually assumes that the density 
is expanded as 

oo 

pch (r) = Y^ ajn (r) , (24) 

n=l 
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where the basis functions fni'f') are drawn from any convenient complete set and the ex- 
pansion coefficients a^ are adjusted to reproduce the differential elastic cross section. The 
corresponding Fourier transform then takes the form 

oo 

p(g)=i^ch(g) = 5ZaJ„(g), (25) 

n=l 

where 

fn (q) = 47r / dr r^jo (qr) /„ (r) . (26) 

represents basis functions in momentum space. 

Evidently the sum in eq. [21] has to be truncated and this produces an error in the 
determination of the charge density distribution. Another problem is that, as shown by 
eq. [231 the solution of the inverse scattering problem requires an accurate determination 
of the cross section up to large momentum transfers. Electron scattering experiments in 
electron-ion colliders will be performed within a limited range of q and this will produce an 
uncertainty in the determination of the charge density distribution. As we have discussed 
in previous sections, the measurements have to encompass the first few minima of the cross 
section for heavy nuclei in order that the density profile can be mapped. 



We consider two bases that have been found useful [37[ in the analysis of electron or proton 
scattering data. The present discussion is limited to spherical nuclei, but generalizations to 
deformed nuclei can be done. The Fourier-Bessel (FB) expansion (i.e. with /„ taken as 
spherical Bessel functions) is useful because of the orthogonality relation between spherical 
Bessel functions 

rRmULX 1 

/ dr r'^ji (qnr) ji (qmr) = --R^ [jl+l (^n-Rmax)]^ f^nm, (27) 

Jo ^ 

where the g„ are defined such as 

31 (qnRm..) = 0. (28) 

The FB basis implies that the charge density Pch{r) should be zero for values of r larger 
than -Rmax- For example, the basis can be defined as follows 

47r (-l)"i?max 
fn (r) = jo {qnr) 6 (i^^ax -r), /„ (g) = ^ T^^^JO (<?^max) , (29) 

where 9 is the step function, -Rmax is the expansion radius and g„ = nn/Rmax- 
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In principle it is possible to obtain the an coefficients measuring directly the cross section 
at the Qn momentum transfer. If the form factor ([2]) is known at g„, the coefficients a„ can 
be obtained inserting (p7|) and (12^ in the definition ([2]) of the form factor, leading to 

a„ = 2- (.30) 

In general the cross sections are measured at q values different from g^. Using the 
expansion fl29p of the charge density one finds for the form factor the relation 

A ( T\^ 

Fch (q) = — y~] CLn^ 5" sin (gi?max) • (31) 

By fitting the experimental Fch{q) one obtains the a„ parameters and reconstruct the nuclear 
charge density. Not all a„'s are needed. Since the integral of the density, or F(g = 0), is 
fixed to the charge number there is one less degree of freedom. Also, densities tend to zero 
at large r. Thus another condition can be used, e.g. that the derivative of the density is 
zero at -Rmax- Thus, when we talk about n expansion coefficients one means in fact that 
only n — 2 coefficients need to be used in eq. [311 For experiments performed up to g^ax the 
number of expansion coefficients needed for the fit is determined by n^ax — ?max-Rmax/7r. 

A disadvantage of the FB expansion is that a relatively large number of terms is often 
needed to accurately represent a typical confined density, e.g. for light nuclei. One can use 
other expansion functions which are invoke less number of expansion parameters, e.g. the 
Laguerre- Gauss (LG) expansion, 

Ur) = e-'L]!' {2a') , and Uq) = 4^/'(3' {-iT e'^' L^' {2^') , 

where a = r//5, 7 = q(3/2, and L^ is the generalized Laguerre polynomial. Another possibil- 
ity is to use an expansion on Hermite (H) polynomials. In both cases, the number of terms 
needed to provide a reasonable approximation to the density can be minimized by choosing 
/3 in accordance with the natural radial scale. For light nuclei /? = 1 — 2 f m can be chosen, 
consistent with the parametrization of their densities. Then the magnitude of a„ decreases 
rapidly with n, but the quality of the fit and the shape of the density are actually indepen- 
dent of (3 over a wide range. As shown by application to a few cases, the main information 
on skin and halo sizes can be obtained using the FB expansion without problems. 

Figure [7] shows the charge form factor squared for a Fermi function distribution with 
half-density radius R = 5.2 fm and diffuseness a = 1.2 fm (solid line). The dotted fine is 
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FIG. 7: Form factor squared for a Fermi function charge distribution with half-density radius 
R = 5.2 fm and diffuseness a = 1.2 fm (solid line). The dotted line is obtained with the FB 
expansion, eq. [311 up to n = 6. The dashed curve uses up to n = 8. 

obtained with the FB expansion, eq. [311 up to ra = 6. The dashed curve uses up to n = 8. 
One clearly sees that the latter improves the fit to the form factor up to the third minimum 
by increasing n = 6 to n = 8. We use -Rmax = 15 fm, so that adding the n = 8 term improves 
the fit of the distribution including the peak at g ~ 1.75 fm^^, as seen in the figure. 
For real data, the expansion coefficients a„ are obtained by minimizing 



X 



E 



Vi -y{qi, «« 



CTi 



where y{qi, an) is the fitted value of the cross section (form factor) with a set of coefficients 
a„ and yi are the experimental data at momentum g^ with uncertainty (Tj. 

To check the limitations of this procedure, we generate a set of pseudodata for ^B. The 
calculated charge form factor of figures [H and [6] is used to generate 40 data points equally 
spaced by Aq = 0.02 fm~^. These data were given an uncertainty ai linearly increasing 
with Qi, from 1% for g^ = 0.02 fm^^ to 20% for qi = 0.82 fm"^ 

The best fit, with -Rmax = 10 fm, is obtained with n = 10 expansion coefficients. Increas- 
ing the number of coefficients does not improve the quality of the fit, as is shown in figure [8] 
for n = 14. It only produces more oscillations of the density. The reason is that terms with 
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FIG. 8: Hartree-Fock charge density of '^B (solid line). The dashed and dotted lines are the 
solution of the inverse scattering problem using a Fourier-Bessel expansion with n = 10 and n = 14 
terms, respectively. 

larger n's are only needed to reproduce the data at larger values of momentum transfer, as 
shown in fig. [71 The fit to the data for q < gmax is not affected but the presence of these 
new terms introduces oscillations in the charge distribution. A possible fix to this problem 
is to include pseudodata in addition to experimental data. This method is well known in 
the literature 37|]. The pseudodata are used to enforce constraints and to estimate the 
incompleteness error associated with the limitation of experimental data to a finite range of 
momentum transfer. 



V. CONCLUSIONS 

In this work I have studied the electron scattering off light unstable nuclei. This work 
is complementary to previous works in this area (see e.g. ref. 10|, lul, 112, ll3|). Particular 
attention was given to the effect of the neutron (proton) skin on the scattering form factors. 
It was shown that the position of the first minimum is arguably the best signature to look 
for noticeable changes in the charge radius size. 

The evidence of a proton halo is not so clear as in the case of other probes. For example, 
it is well known that ^Be fragments arising from the proton knockout of ^B projectiles 
display a distinctively narrow momentum distribution characteristic of a long tail of the 
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valence proton in ^B 2J]. This feature is a consequence of the peripheral character of the 
knockout process, which is ideal to probe the tail of the bound states. Nonetheless, due to 
the Coulomb barrier the bulk of the charge in the nucleus is confined close to its center. 
Electron scattering is sensitive to this bulk charge and therefore does not display such a 
strong halo signature. 

The minimum information obtained with electron scattering in electron-ion colliders will 
be the rms charge radius, (r^^) . This information "per se" is very valuable. It is sensitive 
to the skin size in a heavy nucleus. But it also depends on the accuracy with which this 
quantity can be measured. As we have shown in the previous sections, it will be necessary 
to go beyond the first minimum to extract information about the central value of the charge 
density. 

Accurate measurements at large momentum transfers are crucial if one wants to describe 
the matter distribution with confidence and have a good comparison with predictions of 
different theoretical models. I have shown with a few basic examples using the Fourier- Bessel 
expansion method that, whereas the matter distribution within the halo is well probed by 
measurements at small momentum transfers, the details of the central distribution requires 
measurements at large g's where inelastic processes may play an important role. This again 
imposes constraints on the information that can be extracted from elastic electron scattering 
off halo nuclei. 

VI. APPENDIX 1 - ANALYTICAL FORM FACTORS 

Here I summarize the analytical expressions for the charge form factors which are used 
in this work. For the nuclear charge density given by a hard sphere (uniform distribution 
with sharp cutoff at R) 

p{r) = {'' - (32) 

0, otherwise, 



the form factor is 



F{q) = -^ [sin (qR) - qR cos (qR)] . (33) 

r/ai 



For an Yukawa distribution, p (r) = po e ^'""{a/r), one gets 



4:71 poO^ 
1 + a2g2 



^(^) . . .2.2- 
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The symmetrized Fermi distribution 39 1 



p{r) 



Po 



cosh{R/a) 



cosh.{R/a) + cosh(r/a) 



leads to the form factor 401]: 

Air'^PQa cosh.{R/a) 



F{q) 



Rcos{qR) lid sm{q R) cosh^n qa) 



q sinh(i?/a) [sinh(7rga) sinh^ (irqa) 

The above expression is composed of oscillating terms damped by exponentials. This is 
better seen taking the limit for qa ^ 1: 

iir'^Poa 



F{q) 
The Fermi distribution 



Q 



[Rcos{qR)-7rasm{qR)] e"^'"'. 



Pir) = Po/ [1 + exp {{r- R) /a}] 



(34) 



(35) 



with central density po, radius R, and diffusiveness a, gives a good description of the 
densities of heavy nuclei. This distribution can be fitted by the convolution p (r) = 
J d^r'pp (r) Pa (r — r') of a hard sphere for pA (r — r') and an Yukawa function 38|] for pp (r). 
The advantage is that the form factor factorizes, F{q) = Fp{q)FA{q), and can be calculated 
analytically as 



F{q) 



Anpo 



sin (qR) — qR cos (qR)] 



1 + q^al 



(36) 



The term inside the first parenthesis comes from the hard sphere (uniform distribution 
with sharp cutoff at R) and the second parenthesis is the damping term due to an Yukawa 
diffuseness surface with width ay- Figure [9] compares F{q) of eq. [36]with that obtained with 
the numerical integration of a two-point Fermi density distribution, eq. [351 The parameters 
for Al {R = 3.07 f m , a = 0.519 fm), Cu {R = 4.163 f m , a = 0.606 fm), Sn {R = 5.412 
f m , a = 0.560 fm), Au {R = 6.43 fm , a = 0. 541 fm) , and Pb {R = 6.62 f m , a = 0.546 
fm) (to simplify the figure, 1 did not plot the curves for Au and Cu). In all cases I used for 
the Yukawa function parameter in eq. (136|) ay = 0.7 fm. We see that the agreement is very 
good. 

For light nuclei, it is more appropriated to use Gaussian densities, with the Gaussian 
parameter a = w | (r^)^^, where \/Jr^y^ is the root-mean-square radius of the matter 
density. For example, for carbon, a = 2.018 fm. For a Gaussian density parametrized as 
p (r) = Po exp (— r^/a^), one gets 

F(g) = {na^fpoexp {-q^a'/A) . (37) 
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FIG. 9: Comparison between form factors of a convoluted Fermi- Yukawa distribution calculated 
numerically using eq. [35] (solid lines), or with the analytical equation [36] (dashed lines). 



To simulate nodes in the wavefunctions of light nuclei, a sum of gaussian distributions can 
be used, including terms proportional to r"exp (— r^/a^). The form factors arising from 
these terms can be obtained from nth-order derivatives of eq. [37] with respect to 1/a^. 
For halo nuclei the following parametrization can be adopted 



p (r) = pi exp {-r^/al) + pa y—J exp {-r/a2) , 



(38) 



where the first term describes the density of the core and the second term describes the 
extended halo density. The second term blows up as r -^ and has only a meaning in 
the description of the long tail characterizing the halo wavefunction. If the core contains 
Al nucleons and the halo contains A2 nucleons, the form factor for the above distribution 
becomes 



A, 



F{q) = Al exp i-q'al/A) + -— ^ 



+ a^q 



2' 



(39) 
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